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1. Introduction 

The general relativity as the theory of gravitation is triumphantly tested in "classical 
experiments" [Ij on 

• the deflection of light by the Sun, 

• the perihelion precession of Mercury, 

• the gravitational redshift of electromagnetic radiation, 

• the time delay of signal from satellites due to the curved space-time, and 

• the gyroscope precession during the orbital motion around the Earth as caused by 
the spatial curvature [21 [3]. 

In addition, an expanding Universe is the prediction being inherent for the general 
relativity. Sure, the general relativity conceptually is the perfect theory of classic gravity. 
In this respect, we usually expect that it is valid, indeed, until effects of quantum gravity 
would be essential at Planckian scales of energy that is unreachable in practice. However, 
such the point of view, perhaps, is actually broken: 

The general relativity itself gives us a brilliant tool in order to search for indications, 
which signalize on breaking down its validity: while observing a motion by inertia, we 
get a curvature of space-time, which can be inserted into Einstein equations, that yields 
a tensor of energy-momentum for an appropriate substance, and if properties of the 
substance are mysterious and unpredictable, then we get a hint for suspecting of the 
incorrect description for the nature. This is exactly the case of hypothetic dark matter 
(see, for instance, review in f5]): in the framework of general relativity, it should be 
inevitably introduced as a transparent pressureless substance dynamically isolated from 
the ordinary visible matter made of known, well studied particles, except the interaction 
via the gravity, so that properties of dark matter are artificially tuned. This tuning has 
various aspects. 

First, rotational curves in disc-like galaxies, i.e. dependencies of rotation velocities 
of stars versus a distance to the galaxy center, if described by the law of Newtonian 
gravity, requires the introduction of dark matter with a tuned spatial distribution. 
Unexpectedly, the dark matter halo is inessential in regions, wherein the gravitational 
acceleration caused by the visible matter, is greater than a critical value while the 
halo starts to dominate in regions, wherein the acceleration by the visible matter is less 
than cjQ [H El El [7] . That was M. Milgrom who first introduced the critical acceleration 
in the description of rotational curves [8] . It is spectacular that the critical acceleration 
is universal: it does not depend on a genesis of disk-like galaxy, and it is the same 
for any studied disc-like galaxy. Unbelievably, an amount and spatial distribution of 
dark matter is tuned to the amount and distribution of visible matter in order to 
form in cosmic collisions the dark halos in disc-like galaxies with the same universal 
critical acceleration. In the framework of general relativity, there is no straightforward 
dynamical mechanism for a deduction of such the universal acceleration. Anyway, the 
deduction looks to be very artificial, most probably, it certainly could be the fine tuning. 
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The critical acceleration subtly binds the dark matter to the visible matter. If this 
relation is dynamical, then it is not due to Newtonian gravity, that does not include 
any critical acceleration. Further, in deep regions of dark halo dominance, the rotation 
velocities tend to constant values fo, that empirically satisfy the baryonic TuUy-Fisher 
law [9]: 



where M stands for the mass of visible matter in disc-like galaxy, G denotes the 
Newton constant. Again, the dark matter halo is tuned, so that the star motion 
within the halo strictly correlates with the usual visible matter, while the constant 
is universal [10]. Finally, in order to complete the first item of argumentation, features 
in distributions of visible matter, no doubt, are imprinted in rotational curves even 
in regions of dark matter dominance hence, features of dark matter distributions 
are tuned to the visible matter, though we have no dynamical reasons for such the 
correlations in the framework of general relativity. Moreover, the morphology of spatial 
distributions is absolutely different for the baryonic and dark matter in disc-like galaxies: 
an exponentially falling central bulge and thin disc of stars and gas in contrast to 
power-law decline of dark spherical halo. Thus, the universal critical acceleration is the 
mysterious quantity for the general relativity, that cannot be predicted, while its notion 
emerged empirically. The critical acceleration reveals the fine tuning of hypothetic 
dark matter to the ordinary visible matter. 

Second, in cosmology with the observed accelerated expansion of Universe by 
data on a dependence of brightness of type la supernovae versus the red shift 
[TT| [T2| [T3t [TU [15] , the general relativity was have to introduce the extended dark sector, 
which includes a homogeneous dark energy in addition to the nonhomogeneous dark 
matter. In the simplest case, the dark energy can be represented by the cosmological 
constant, otherwise it should be described by a homogeneous fluid X with the state 
parameter wx being the ratio of pressure px to energy density px, wx = Vx/ Px close 
to vacuum value of —1, in contrast to the pressureless dark matter with wdm = 0. 
Evidently, the nature of dark matter and dark energy is very different. But surprisingly, 
the energy density of dark energy, or the value of cosmological constant A, is finely 
tuned to the critical acceleration [1], so that 



Therefore, the dark energy should be inherently connected to the unexpectedly 
correlated dynamics of dark matter and ordinary matter. But the coincidence of ([2]) is 
mysterious for the general relativity. 

Nevertheless, the general relativity applied to the cosmos still looks formally viable 
in the form of concordance model: the ordinary matter balanced with the cold dark 
matter (CDM) and cosmological constant A in the fiat space, the ACDM variant with 
a spatial curvature compatible with zero in limits of uncertainties. Moreover, there are 
two important successes in the model: 1) a correct fitting for observed anisotropy of 
cosmic microwave background radiation (CMBR) [16] that becomes possible due to a 



= GM~go, 



(1) 



Gpx ~ A ~ 



(2) 
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tuned amount of non-baryonic dark matter, and 2) an appropriate baryon to photon 
ratio consistent with a current status of big bang nucleosynthesis (BBN) (see review by 
Fields B D, Sarkar S in [H]). 

Indeed, the anisotropy of CMBR is suppressed by 5 orders of magnitude with 
respect to the CMBR temperature, and it is caused by a propagation of sound waves in 
a hot photon-electron-baryon medium up to a moment, when the electrons bind to nuclei 
to form neutral transparent gas. The snapshot of Universe at the time of decoupling 
of photons evolves to us, and it represents acoustic peaks in the following multipole 
spectrum of temperature fluctuations 



where ni^2 denote directions in the celestial sphere. Pi are Legendre polynomials of 
multipole number /. The spectrum, i.e. C/, depends on 

• the Universe evolution, 

• a primary spectrum of inhomogeneity, and 

• the propagation of inhomogeneity during the evolution. 

Then, the Universe evolution is well described by ACMD |16] and it can be extrapolated 
to the age of Universe, when the snapshot of CMBR was done, i.e about 380 thousands 
years after the big bang. The primary spectrum of inhomogeneity is suggested to be 
close to the Harison-Zeldovich distribution of so called "no-scale" limit at a spectral 
index ns{k) = 1: 



where Sp{r) = p{r) — p denotes a contrast of energy density at comoving coordinate, 
k is a wave vector conjugated to the comoving coordinate, A stands for an amplitude 
at a reference value of k^. The spectral index and amplitude are subjects to fit the 
observed spectrum of Ci. Finally, the propagation of inhomogeneity includes the sound 
and further smearing of waves by the gravitation. The concordance model of cosmology 
in the framework of general relativity well fits Ci with the fiat space and tuned amount of 
dark components [16]: relative fractions of dark matter Qbu ~ 20% and dark energy in 
the approximation of cosmological constant Qa ~ 76%. The baryonic matter composes 
only Qh ~ 4%. This value is dictated by heights of distinct initial three acoustic peaks 
in the multipole spectrum, C/. That is the dark matter fraction, which regulates the 
relative heights and positions of peaks up to small variations due to the parameters of 
primary spectrum of inhomogeneity. 

Next, the amount of baryonic matter and the temperature of CMBR fixes the 
baryon-to-photon ratio of densities r]h = ni,/n^, that is the only free parameter in 
calculation of elements abundance during the big bang nucleosynthesis. The current 
state of measurements of elements abundance extrapolated to the primary abundance 
is compatible with r]h extracted from the ACDM fit of CMBR anisotropy [IT] . 




(3) 




(4) 
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Thus, the success of concordance model stimulates direct searches for an appropriate 
heavy dark matter particle at colliders and underground big-volume experiments, 
sensitive to suppressed, but non-zero cross sections of dark matter interaction with 
the ordinary matter. 

However, even a discovery of candidate for the dark matter particle would not 
withdraw the problem of ad hoc tuning of dark sector. Moreover, it would sharpen the 
need to search for the dynamical reasons causing the adjustment of dark matter, i.e to 
look beyond the general relativity. 

A model of gravity involving the critical acceleration should naturally include 
both empirical laws such as the Tully-Fisher relation at the galactic scales and correct 
descriptions of Universe evolution, observed features of CMBR, large scale structure 
and elements abundance, so that the model would give a successful approach being 
alternative to the general relativity in cosmology, of course. At galactic scales, 
M. Milgrom invented the modified Newtonian dynamics (MOND) [8J stating the 
gravitational acceleration g 

9 Cigf-go) = - V0M, Ciy) = S^ , (5) 

yi + r 

where the critical acceleration extracted from the modern analysis of rotational curves 
is given by = (1-24 ± 0.14) x 10^^° m/s^ [10], while (pM denotes the Newtonian 
gravitational potential of ordinary matter, satisfying V^</)a/ = —An G pm. At g/go ^ 1, 
we get the Newtonian limit of gravitational forc^, while at g / go 1 the Tully-Fisher 
law is satisfied by construction. Note that ([5]) successfully predicts the rotational curves 
by the given distribution of visible mattei[§| with appropriate imprints of its features, see 
review in |4]. 

However, the straightforward insertion of ([5]) into the dynamics at cosmological 
scales would results in an inconsistent distortion of vacuum homogeneity during the 
evolution of Universe, for instance, i.e. in the vacuum instability as was shown in [T8] . 
wherein authors offered to introduce the cosmological behavior of critical acceleration 
in the form of 

^0 ^0 = ^0 la^l, (6) 

where the distance is determined by comoving coordinate r and scale factor of Universe 
expansion a{t), so that x = a{t) r. Then, the homogeneous cosmology with the gravity 
law modified at accelerations below the critical value of Qq is consistent, that constitutes 
the cosmological extrapolation of MOND [l_8j. The cosmological regime is matched 
to MOND at a size of large scale structure \x\\ss, i-e. at the characteristic scale of 

I Sub-leading terms are suppressed, so that the force at the Earth and in the Solar system is not 
distinguishable from the Newtonian one. 

§ The only parameter of fitting the rotational curves within MOND is the light-to-mass ratio, which 
strictly correlates with astrophisical expectations for given galaxies. Moreover, in gas-rich galaxies this 
uncertainty is absent, that means the MOND predicts rotational curves with no adjustment of any 
parameters, see details in [4]. 
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inhomogeneity, that is related to the acoustic scale and sound horizon in the baryon- 
electron-photon plasma (see details in [T8]). 

In the framework of cosmological extrapolation of MOND with the interpolation 
function C(?/) in (|5]) the gravity equations for the evolution of homogeneous and isotropic 
Universe can be written in the form [18] 

{Rli'Q ' = ( ' + ') ' , (7) 

in terms of Ricci tensor R^j^ for the metric rj^y and 4-vector in the direction of 
cosmological time = (1, 0), wherein the matter energy-momentum tensor T^^, defines 

i?;;M = 87rG(T;-ir/;:T), (8) 

while the extra tensor of curvature K'^ is the Ricci tensor of de Sitter space-time being 
both homogeneous in time and space as well as isotropic 

k; = Hv;- (9) 

In addition to ([7]) the conservation of energy-momentum V ^Tjf = is hold, of course. 

Then, in the limit of general relativity we put {R'j^^^C.u) ^ {^ll^^^u) , and we find 
R'^ ^ R'^, that results in the Einstein equations 

R-; = 8nG [t; - ir^;:T) . (10) 

In the limit of {R^^^^^v) ^ {^^^'^^u) we get the modified evolution of Universe, 
effective at present. 

It is important that parameterizing the size of large scale structure by |£c|iss ~ A^/i^o 
at ~ Hq/X and moderate value of A ~ ^, we find the simplest solution for the 
coincidence problem, because the Milgrom acceleration c/q ~ XHq becomes close to the 
scale of cosmological constant A ~ Hq. Moreover, we will see that the dark matter 
fraction of energy is also regulated by the value of A. 

Eq. ([7]) successfully fits the evolution of Universe measured by observing the 
brightness of type la supernovae versus the red shift [18]. So, the stellar magnitude 



fl = fiabs + 5\og-i^Q dL{z) +25, (11) 

depends on the photometric distance di (in Mpc), determined by the Hubble constant 
evolution H{z) = a/ a, 

z 

f dz 

dL{z) = {l + z) j (12) 



where iiabs is an absolute stellar magnitude of light source at the distance of 10 pc. We 
show the Hubble diagram for the type la Supernovae in Fig. [TJ The mean deviation 
squared per degree of freedom gives x^/d.o.i. = 1.03 for our fit with the following 
assignment of parameters: 

go = -0.853, zt = 0.375, 

h = 0.71, fib = 0.115, ^ ' 
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Figure 1. The magnitude of supernova luminosity versus the redshift z. The data 
with error bars are taken from the Union2 collection [T5] and the results of Hubble 
Space Telescope on the type la supernovae [12] ■ The curve represents our fit in the 
framework of modified cosmology. 



where go determines the deceleration parameter at red sift z = 0, q{z) = —d/aH^{z), 
Zt stands for the red shift, when the acceleration is equal to zero, h parameterizes the 
Hubble rate at z = via Hq = h ■ lOOkm/s ■ Mpc~^. Accepting the prescription of 



9'o- 

we can find that 
go - 



KoH, 



05 



i^O^^, ((1 + Ztf - 1) 



(14) 



(15) 



when the energy density determined by the cosmological constant, is given by the energy 
budget of Universe, f^A = 1 — Qb- The same values of deceleration, go, and red shift of 
transition from the deceleration to acceleration, Zf, could be obtained in concordance 
model of ACD^iJ at 

= (16) 

lib go 

Therefore, the modified gravity with the critical acceleration, i.e. at nonzero Kq ~ 7.9, 
determines the ratio of matter to baryon fractions of energy, so that at go ~ —1 we get 
QM/^b ~ Kq, which correlates with the scale of large scale structure, considered above, 
Ko ~ 1/A. 



II The ACDM parameters are marked by bars. 
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Note that the quahty of fit is not sensitive to valuable variations of baryonic fraction 
while it was extracted from the appropriate value of sound horizon 

t{,z) 

rs{z) = j Csdt, (17) 



where the speed of sound is determined by the baryon-to-photon ratio of energy i?, so 
that 

The value of sound horizon at 2; = 0.2 and z = 0.35 was extracted from baryonic 
acoustic oscillations (BAO) |19], which, in the case of baryon matter only, favor for the 
enhanced estimate of Qb shown above. The same conclusion follows from the calculation 
of acoustic scale in the spectrum of CMBR anisotropy [16], 



(1 + z^)rs{z^) 

here z^, is the redshift of decoupling, when due to the recombination of electrons with 
protons the medium becomes transparent for photons (see analytical approximations for 
in terms of baryonic density, matter density and Hubble constant in [20].). WMAP 
gives Ia = 302.69±0.76, while we deduce Ia = 302.5 [18] compatible with the uncertainty 
of measurement. Hence, we expect that scale features of CBMR anisotropy could be 
fitted in the framework of cosmological extrapolation of MOND. 

In present paper we describe the procedure of fitting the CMBR anisotropy 
spectrum with the model of modified gravity and point to accepted approximations 
in Section [2l Some actual problems associated with the theory and phenomenology of 
our model of modified gravity, are considered in Section |3l We present an itemized 
discussion of model verification and falsification in Section |H 

2. Fitting the CMBR anisotropy 

The tool for the calculation of multipole spectrum of CMBR anisotropy [211 1221 EH] 
operates with the Friedmann equation, which is not valid in the framework of 
cosmological extrapolation of MOND. Nevertheless, we can integrate out the dynamical 
equations of our model in order to extract the Hubble rate at any red shift and to 
parameterize it with a mimic dark energy contribution in addition to the fraction of 
baryonic matter. Indeed, exploring the general relativity in the case of baryonic matter 
and dark energy in the isotropic homogeneous curved space we get 

^''/''-^'+^^x(a)(l + 3^x(a)), (21) 



Hi a3 
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where Vt^ stands for the contribution of space curvature, and the energy budget holds 
VL}j + VtxiX) + fifc = 1- Excluding VLx{.ci), we get the expression for the dark equation of 
state 



Here we insert the expression for the acceleration that follows from eq. ([7]) for the 
modified gravity including the actual values of energy fractions for baryons and vacuum, 
Vth and f^A, respectively. 



Note that the red shift of transition from the acceleration to deceleration of Universe, 
Zt is related to the fraction of cosmological constant due to 



Again, eq. fl23|) clearly shows that putting g'^ equal to zero, we get the dynamics of general 
relativity, otherwise near a = the dynamics enters the region of strong dominance of 
modification and the greatest deviation from the general relativity, as it is actual at 
present. 

Integrating out fl23|) at the initial condition a/a{t = tg) = Hq, we obtain0the scaling 
quantity H/Hq required for the complete definition of r.h.s. in ( !22|) . Fig. |2] shows the 
behavior of wx versus the scale factor a as we have calculated in the cosmological 
extrapolation of MOND with parameters listed in ( IT3|) at = 0. Small variations of 
the spatial curvature in limits \flk\ < 0.02, deceleration parameter qo, transitional red 
shift Zt and baryonic fraction Qi, within 10% lead to negligible changes, which are only 
just visible in the figure. We emphasize that the modification of gravity predicts the 
very specific dependence of equation of state for the dark energy, that we will discuss 
in Section m 

After the definition of referenced homogeneous evolution of Universe under the 
modified gravity, we can use the standard tool for the calculation of CMBR anisotropy 
spectrum [211 |22l [23]. However, in this way the propagation and smearing of sound 
waves would be described in the framework of general relativity with no dark matter, 
i.e. at Qbm = 0, while we have to modify this procedure in accordance with a structure 
formation under the modified law at accelerations below the critical one. The problem 
is that such the modification of perturbation transfer function is not linear, and the 
appropriate machinery of calculation is not yet developed. That is missing point of our 
consideration, of course. Nevertheless, the extensive usages of tool have shown that the 
influence of dark matter on the spectrum is reduced to relative enhancement of third 




(22) 




(23) 




% In practice, we use the scaling variable r — tHo, that completely covers the differential equations 
under consideration. 
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peak, whereas this enhancement is due to enforcing the gravity. So, since the modified 
gravity produces the very similar effect of enforcing the gravity, we can expect that it 
could results in the same fine feature as concerns for the enlarging the third peak. 

In this respect, we have to emphasize that a complete axiomatic approach based 
on a formulation of action for a modified gravity in terms of given extended set of 
gravitational fields has got the advantage in calculating of various predictions including 
the CMBR anisotropy. So, we can mention the following fully relativistic schemes (more 
examples and references find in [4j): 

• Bekenstein's theory of tensor-vector-scalar (TeVeS) gravitational fields j24] 
involving Maxwellian vector field and reproducing the critical acceleration in the 
case of isolated gravitational source, equivalent to MOND, 

• Moffat's modified gravity [251 EE! [27], giving the approximation of gravity law 
similar to the MOND, 

• generalized TeVeS theories [28] with non-Maxwellian vector field. 

However, first, these theories include dark gravitational fields actually replacing the 
dark matter that looks like a refinement of problem. Second, all of them have the strict 
theoretical illness: there are configurations with unlimited, infinite negative energy, that 
leads to instability of physically sensible solutions (see details and references in [1]). 
Third, the critical acceleration is still introduced ad hoc with no reasonable connection 
to the present Hubble rate or the cosmological constant. 

So, we prefer for the phenomenological approach, which does not introduce new 
artificial and heuristic notions. In this way, we can investigate the role of critical 
acceleration in the modified gravity by studying various phenomena step by step in 
order to find fundamental features and differences from the general relativity. 

Wx 

0.2 

0.0 
-0.2 
-0.4 
-0.6 
-0.8 
-1.0 
-1.2 




Figure 2. The equation of state wx versus tiie scale factor a in tlie model of modified 
gravity involving the critical acceleration. 
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At present, we try to fit the spectrum of CMBR anisotropy by using the 
modified evolution of homogeneous Universe and optimizing the primary spectrum of 
inhomogeneity, which develops as the sound smeared by ordinary gravity. So, we suggest 
that a modification of smearing will not be very crucial for the main features of spectrum. 

The results of such the fitting are presented in Figs. [3H5l First, we study the fit of 
WMAP data [16] with the Harrison-Zeldovich primary spectrum (HS) at = 1, which 
is shown in Fig. |3] by dashed line. It is spectacular that the no-scale HS prescription 
correctly reproduces the angular scale of multipole momentum, i.e. the position and 
profile of first acoustic peak. This feature is obtained due to the correct adjustment of 
this scale by the sound horizon in the case of no dark matter [18] as mentioned in the 
Introduction. Then, we find that the running of spectral index ns{k) = n^^^ + n'^ In fc/Zco 
leads to suitable description of both, first and second acoustic peaks in the modified 




Figure 3. The spectrum of CMBR anisotropy calculated in the model of modified 
gravity, l{l + l)C;/27r (in ixK^) versus the multipole number I, in comaprison with 
WMAP data [16]. The dashed line represents the Harrison-Zeldovich approximation, 
while the solid curve corresponds to the running spectral index (as described in the 
text). 
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cosmology with no dark matter. Setting parameters equal to the following value^: 





= 1.625, 


K 


= 0.24, 




A 


= 5.5 X 10"^ 




= -0.015, 


(24) 


Zre 


= 23, 


ko 


= 0.04Mpc"\ 





we get the fit shown by solid line in Fig. |3l Set f l24|) needs comments. 

First, the running of spectral index allows to adjust the relative height of second 
acoustic peak in the spectrum0- This running is dynamically essential and numerically 
significant. Moreover, because of sizable value of slope of spectral index, we expect that 
the approximation neglecting the higher orders of expansion versus the logarithm of 
comoving wave vector Ink, would be inaccurate at large intervals of multipole moment. 

Second, we introduce a small spatial curvature in order to adjust the position of 
peaks, which have been displaced under the strong running of spectral index. Then, the 
spatial curvature essentially improves the quality of fit. Remember, that such the low 
value of spatial curvature is consistent with the equation of state we have deduce from 
the modified gravit}!^. Moreover, acceleration a does not involve the spatial curvature 
because of its specific value of state parameter Wk = — |, giving pk + 3pk = 0. We 
note also that a nonzero value of spatial curavture obeys the scaling If^Q^fel ~ 1, that 
is a feature consistent with our solution or treatment of coincidence problem (see the 
Introduction). 

Third, the red shift of reionization Zre, when inhomogeneities of cold g clS SbTQ 
contacted due to the gravity in order to form stars, refers to the heating of gas, that 
transforms it to plasma again. It causes a reduction of intensity of radiation passing 
through the hot secondary plasma. So, the amplitude of primary CMBR, A correlates 
with the red shift of reionization. In MOND with the enhanced gravity at the stage of 
diluted gas, one expects that the star formation starts early than in the general relati vity 
[4j, hence, we try to improve the fit quality by enlarging the red shift of reionizatio nffl . 
The corresponding improvement of mean deviation squared is equal to -0.11 per degree 
of freedom. The optical depth is shifted from 0.9 to 0.78. So, we refer this adjustment 
as fine effect beyond a sole significance, as well as slow variation of pivot wave vector 
ko from 0.05 to 0.04 Mpc~\ 

Thus, in the simplest way of modifying the background cosmology by the 
cosmological extrapolation of MOND, we find that WMAP data [16] can be fitted at 
the following mean deviation squared per degree of freedom: 

xVd.o.f. = 1.34, (25) 

+ We list the quantities, which we change from default values assigned in the tool [1T1[331[53]. Other 
quantities have been set to its standard prescriptions. 

* The opportunity to fit the second peak in the model with suppressed dark matter was considered in 

m- 

H Note that at a the equation of state for the mimic dark energy tends to wx{0) ~ — 5, which 
points to a possibility of nonzero spatial curvature. 

ft We set Zre via "a sizable change": the typical value of z^e = 11 in ACDM has been enlarged twice 
for a distinguishability with no strict reasons or prerequisites. 
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Figure 5. The same as in Fig. 21 but focused at the region of third and further peaks 
{I > 600). 
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which is just 0.19 worse than the typical value of x^/d-o.f. in the concordance model of 
ACDM. Fig. m clearly shows that at / > 600 the WMAP data suffer from huge statistical 
and systematical uncertainties. Indeed, the inclusion or exclusion of data at / > 600, 
in practice, do not change x^/d.o.f. Therefore, the WMAP data beyond the first and 
second peaks are not conclusive, at all to the moment. 

Then, higher multipoles can be studied due to the BOOMERANG and ACBAR 
data sets [301 [31], [32] . We present the comparison in Fig. [51 Evidently, the dark matter 
indication is reduced to the different form of third peak rise. Therefore, we can hope 
that a nonlinear smearing of sound waves in MOND could give the same effect in the 
region of third peak, too, although to the moment we see the clear tension between the 
data in the region of third acoustic peak and the simplest version of fit based on the 
modified cosmology of homogeneous Universe only. Thus, we need the development of 
procedure to calculate the transfer function of inhomogeneity versus the red shift [33] 
in the case of nonlinear MOND. The same is true as concerns for the simulation of large 
scale structure formation visible at present time. 

Nevertheless, we insist on quite the successful fitting of CMBR anisotropy in the 
framework of cosmological extrapolation of MOND, as the variant of modified gravity 
involving the critical acceleration in the case of homogeneous matter. 

3. Main problems 

The enhanced value of baryon fraction in the energy budget, ^ 0.115 implies the 
enhanced value of baryon-to-photon ratio 77?, being the only free parameter of BBN 
calculation [17j. The BBN takes place during the period when the deviations from the 
general relativity are negligible. Therefore, the primary abundance of light elements in 
cosmological extrapolation of MOND should differ from the BBN estimated within the 
concordance model of ACDM. 

At present the data on the helium abundance has large uncertainties which are 
compatible with both models under consideration. Then, the deuterium and lithium 
primary abundances are clearly able to discriminate between the enhanced value of rj^ 
and its concordance value. However, these quantities are not measured directly, they 
are extrapolations from a suggested model of evolution. 

Indeed, we observe the visible sources of light, that mean the primary matter is 
contracted in stars with further development of nuclear reactions, and a model of star 
evolution has to reproduce the primary values of abundances. In this respect, one 
has to take into account the different red shifts of star formation, i.e. different ages 
of luminous objects, as considered in the framework of general relativity or MOND, 
of course. Next, at present the standard model for the extrapolation to primary 
abundances is not self consistent, because it predicts different ratios of both deuterium 
to lithium-7 and lithium-7 to lithium-6, whereas the former ratio is in bright tension of 
prediction with the extracted values, while the later is in a deep contradiction (about 
three order of magnitude!). Thus, the present BBN status can not be surely conclusive. 
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It is desirable to get more better reliability of procedure for the empirical extraction of 
primary abundance of light elements. 

On the other hand, the doubling of baryonic fraction should appear in cosmological 
effects. So, this doubling could appropriately explain a missing mass in galactic clusters, 
as found even within MOND. Then, a significant portion of baryonic matter should be 
in cold form (Jupiter-like objects), for instance. Note that even in ACDM the visible 
matter composes the tenth fraction of all baryonic matter, only, hence, the most of 
baryons are in cold form. 

Next, the visible large scale structure should be explained by appropriate 
propagation of primary spatial inhomogeneities. However, this issue of modified gravity 
involving the critical acceleration is not yet developed because of nonlinearity of the 
problem. 

4. Discussion and conclusion 

In this paper, we have shown that the cosmological extrapolation of MOND as the 
modified gravity involving the critical acceleration, can successfully reproduce main 
features of multipole spectrum of CMBR anisotropy. 
Let us list conclusions of our investigation. 

(i) The modified dynamics adjusted to empirical values of the Hubble rate, sound 
horizon in the baryon-photon medium, acoustic scale in the multipole spectrum 
of CMBR anisotropy and magnitudes of type la supernovae at red shifts 2; < 2, 
mimics the dark energy with the very specific equation of state wx, shown in Fig. 
[2j This is the falsifiahle prediction of cosmological extrapolation of MOND. It 
can be verified in the nearest future by extensive measuring of type la supernova 
magnitudes versus the red shift. Even z < 2, i.e. the scale factor variation 
within the interval 0.3 < a < 1, would be enough in order to make decision on 
the direct falsification of cosmological model involving the critical acceleration. 
Moreover, such the exotic behavior of dark energy state parameter wx, if would be 
confirmed, will be marginally artificial for the general relativity, that would mean 
the straightforward indication of its inadequateness. 

(ii) The spectral index of primary spatial inhomogeneity ns{k) = n^^^ + n'^lnk/ko, 
essentially runs. 

• The character of running is model- dependent, and 

• It is very different in the concordance model of general relativity and in the 
cosmological extrapolation of MOND: in the modified gravity the running is 
rather fast, and its parameters signalize on the hybrid (multified) inflation, 
that could generate such the spectrum, while the general relativity gives the 
slow running, which preferably corresponds to the simplest top-hill inflation 
due to a single infiaton field [M] . 
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• The modified gravity gives nf^ > 1 at ri^ ~ | > 0, while the general relativity 
results in nf* < 1 at n'^ — )• 0. 

• The fast running probably indicates the need to improve the calculation tool in 
order to include higher derivatives of spectral index with respect to logarithm 
of wave vector. 

(iii) The spectral index of primary inhomogeneity define initial conditions for the 
transfer of inhomogeneity during the evolution, which can be observed in baryonic 
acoustic oscillations in large scale structure of present Universe |T9]. The procedure 
of calculating the transfer function in the framework of modified gravity with the 
critical acceleration is nonlinear, and it is not still developed, that does not allow 
us to make a comparison with data, at present. The structure growth is enhanced 
in MOND, and it can probably need for additional mechanism of smearing the 
acoustic oscillations 

(iv) The modified gravity in our version results in doubly enhanced fraction of baryons, 
approximately. This means that baryon-to-photon ratio is twice large, at least. 

• More reliable estimates of primary abundance of light elements is required 
in order to discriminate the general relativity from the modified gravity by 
the big bang nucleosynthesis. Therefore, BBN can give the falsification of 
cosmological extrapolation of MOND. 

• The doubling of baryons in the form of cold baryonic matter (Jupiter-like 
objects, for instance) should be found in observations. For instance, the 
mass deficit in galaxy clusters described within MOND, can signalize on the 
appropriate enhancement of baryons. 

(v) The multipole spectrum of CMBR anisotropy needs improvements of accuracy in 
the range of third acoustic peak. If the model of modified gravity will still miss 
the correct description of third peak after such the improvement, then this would 
point to the extension of simplest consideration by strict inclusion of inhomogeneity 
propagation within the modified dynamics. 

Finally, we have shown that the coincidence problem of general relativity is 
inherently solved in the framework of cosmological extrapolation of MOND: the critical 
acceleration is connected to the extra Ricci tensor of de Sitter space, involved in 
the gravity equations; then, it is naturally correlates with the cosmological constant. 
In addition, the modified gravity is mostly effective at zero acceleration of Universe 
expansion. That is why the coincidence notion is actual at present. 

Keeping in mind soluble problems mentioned, we state that the cosmological 
extrapolation of MOND is quite successful in cosmology. Moreover, we can falsify it in 
the nearest future, although the same note on the verification is actual for the general 
relativity, too. 
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